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ABSTRACT 

We search for chemically-distinct stellar components in the solar neighbourhood using 
a compilation of published data. Extending earlier work, we show that when the 
abundances of Fe, a elements, and the r-process element Eu are considered together, 
stars separate neatly into two groups that delineate the traditional thin and thick 
disk components of the Milky Way. The group akin to the thin disk is traced by 
stars with [Fe/H] > —0.7 and [a/Fc]< 0.2. The thick disk-like group overlaps the thin 
disk in [Fe/H] but has higher abundances of a elements and Eu. Stars in the range 
— 1.5 <[Fe/H]< —0.7 with low [a/Fe] ratios, however, seem to belong to a separate, 
dynamically-cold, non-rotating component that we associate with tidal debris. The 
kincmatically-hot stellar halo dominates the sample for [Fe/H]< —1.5. These results 
suggest that it may be possible to define the main dynamical components of the solar 
neighbourhood using only their chemistry, an approach with a number of interesting 
consequences. With such definition, the kinematics of thin disk stars is found to be 
independent of metallicity: their average rotation speed remains roughly constant in 
the range —0.7 < [Fe/H] < +0.4, a result that argues against radial migration having 
played a substantial role in the evolution of the thin disk. The velocity dispersion of the 
thin disk is also independent of [Fc/H], implying that the familiar increase in velocity 
dispersion with decreasing metallicity is the result of the increasing prevalence of the 
thick disk at lower metallicities, rather than of the sustained operation of a dynamical 
heating mechanism. The substantial overlap in [Fe/H] and, likely, stellar age, of the 
various components might affect other reported trends in the properties of stars in the 
solar neighbourhood. A purely chemical characterization of these components would 
enable us to scrutinize these trends critically in order to understand which result from 
accretion events and which result from secular changes in the properties of the Galaxy. 

Key words: Galaxy: abundances, Galaxy:kinematics and dynamics, Galaxy: forma- 
tion, Galaxy: evolution, galaxies: structure, galaxics:formation 



1 INTRODUCTION 



ponent truly distinct from the thin disk requires additional 
evidence. 



The discovery that the vertical distribution of stars in the 
solar cylinder is best approximated by two exponential laws 
suggested the presence of two components of different scale- 
heights and ushered in the concept of the Galactic thick disk 
(IGilmore fe Reidll983h . The concept of separate "thin" and 
"thick" disk components in the Milky Way, however, is use- 
ful insofar as they refer to different types of stars; i.e., dis- 
tinct in ways other than kinematics or spatial distribution. 
Since there is no a priori reason why the vertical structure 
of galaxy disks cannot be more complex than a single expo- 
nential, confirming the identity of the thick disk as a com- 



In that regard, the finding that thick disk stars are 
substantially older and more highly enriched in a ele- 
ments than their thin disk counterpar ts supports the two 
component nature of the Galactic disk dGratton et alJ 1996 ; 



Fuhrmannl 19981; jProchaska et al" 200d : Bensbv et al 
vsd 



2003 



Gilmore fe Wvsel Il985l : iReddv et alTl2006l : iBensbv et all 
20071 : lFuhrmannll2008l h The thin and thick disks have metal- 
licity distributions that overlap in [Fe/H] but that differ, at 
given [Fe/H], in their kinematics, age, and a content. This 
dichotomy in properties at fixed metallicity requires adjust- 
ments to traditional models of chemical evolution, which 
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have invoked violent accretion events, as well as episodic 
hiatus in star formation, to explain th e data (see, e.g., 
IChiappini et al.lll997l : iReddv et al.ll2006t ). Such events may 
reset selectively the heavy-element abundance of the ISM, 
disrupting the monotonic trends with [Fe/H] expected in 
simple models of chemical evolution and enabling better fits 
to the data. 

The price to pay is one of increased model complex- 
ity. Did accretion events supply mostly gas or stars to the 
disk? Were stars in the thick disk brought into the Galaxy 
during the event or is the thick disk the stirred-up rem- 
nant of an early thin disk? How many accretion events took 
place? What were the masses of the accreted dwarfs? When 
did these events take place? Accretion models provide only 
rudimentary answers to these key questions and, as a re- 
sult, there is little consensus in the community regarding 
the origin of the thick disk. 

This state of affairs has been highlighted by a num- 
ber of recent papers, which argue that the chemo-dynamical 
evidence for accretion events is weak and that all rele- 
vant data can be explained by reconsidering the impor- 
tance of radial migrati on of both gas a nd stars during the 
evolution of the disk jHaywoodl 120081 ; iRoskar et all 120081 ; 
ISchonrich fc Binnevl l2~009ah. Inspiration for this w ork came 
from the realization by Sellwood fc Binnevl |2002l ) that in- 
homogeneities in the disk do not just perturb disk stars dif- 
fusively, easing them gradually into more and more eccentric 
orbits, but can also transport, resonantly, stars across the 
disk without increasing their eccentricity. Stars on nearly cir- 
cular orbits in the solar neighbourhood could therefore have 
formed at different radii in the Galactic disk. The properties 
of local thin disk stars may thus reflect large-scale gradients 
in the Galaxy rather than different conditions in the so- 
lar neighbourhood at the time of their formation. Using a 
simple but plausible model that includes radial migration, 
ISchonrich fc Binnevl (|2009bh demonstrate that most avail- 
able data for the thin and thick disks can be accommodated 
without invoking a violent accretion event at all. 

Deciding between migration or accretion scenarios re- 
quires identifying patterns in stellar properties that point to 
the presence of truly distinct stellar component and assess- 
ing whether their properties are the result of secular evolu- 
tion or accretion events. In general, secular evolution mecha- 
nisms such as radial migration should lead to increased mix- 
ing, blurring the boundaries between components and leav- 
ing dynamical and chemical imprints different from those 
predicted by accretion scenarios. 

Take, as an example, the wide range in metallicity 
spanned by stars on nearly-circular orbits (the thin disk) 
in the solar neighbourhood. In migration-based scenarios 
the metal-poor tail of the local thin disk is populated by 
stars that formed further out in the Galaxy, while the oppo- 
site holds for stars in the metal-rich tail. Stars at these two 
extremes were born with very different angular momenta, 
and this difference would be preserved in local samples even 
after their orbits have migrated to the solar neighbour- 
hood (|Schonrich fc Binnevl ]2009bl ). Therefore very metal- 
rich (poor) stars would be found in the solar neighbour- 
hood at relatively low (high) rotation velocities. A relatively 
clean prediction of migration models is, then, the presence 
of a negative correlation between mean rotation speed and 
metallicity for stars in the thin disk. 



Such correlations have been difficult to study because 
of the widespread practice of assigning stars to the thick or 
thin disk according to their kinematics. Although this may 
be an expeditious procedure, it imposes an obvious bias that 
precludes searching for correlations of the kind alluded to 
in the preceding paragraph. It would be much more useful 
to devise a purely chemical characterization of the various 
components in the Galaxy. After all, chemistry is a relatively 
stable property of a star that relates it to the conditions 
at the time/place of its birth, whereas its distance to the 
center, vertical height, or spatial velocity are far less durable 
features. 

We explore these ideas here using a compilation of data 
for stars in the solar neighbourhood with good estimates 
of their spatial motions and reliable measurements of their 
heavy-element abundances. We focus on spectroscopic sam- 
ples where the abundance of individual elements, such as Fe, 
Mg, Ti, and Ca, are also available. We use these data to ap- 
portion stars to separate components on the basis of chem- 
istry alone, and use their kinematics to place constraints on 
their possible origins. 



2 DATASET 



We have used the compilation of I Venn et all (12004 ), which 
includes abundance data for solar neighbourhood stars with 
good spatial mo tions, supplemented by data in th e more 
rece nt papers of IReddv et ail (|2006l ). iBensbv et ail (|2005l ) 
andU ISSl'll fc Schuster! (|2010l ). This is a heterogeneous com- 
pilation of data from many different sources. We have gone 
carefully through this compilation in order to remove re- 
peated stars and to bring all velocities to a consistent ref- 
erence frame, taking into account the fact that some papers 
use either right- or left-handed UVW frames. The sample 
spans a wide range in metallicity (—4 <[Fe/H]< +0.5), and 
there are reliable abundance measurements of [Fe/H] and of 
the a elements Mg, Ti, and Ca, for 743 stars. Abundance 
measurements of the r-process element Eu are also available 
for 306 of those stars. We shall us e these two subsamples 
in what follows. As in IVenn et alj (|2004l ). the spatial ve- 
locities assume that the Sun moves in a Galactic reference 
frame with (U,V,W)=(9,232,7) km/s and that the LSR has 
Vlsr = 220 km/s. 



3 RESULTS 

Fig. [T] shows, as a function of the iron content, [Fe/H], an 
index measuring the average abundance of the a elements 
Mg, Ti, and Ca, plus the r-process element Eu for all stars 
in our sample. The index is constructed by simply averaging 
the Mg, Ti, Ca, and Eu abundances. As usual, all quantities 
are shown in logio units normalized to solar values. Like the 
a elements, Eu is thought to form mainly in massive stars, 
and it therefore may be combined with the a elements to 
define an index that gauges the importance of massive-star 
nucleosynthesis in a given stellar component. 
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Figure 1. Iron abundance vs an index constructed by averag- 
ing the abundance of the a elements (Mg, Ti, and Ca) plus the 
r-process element Eu. Note the presence of two well-defined "fam- 
ilies" of stars each spanning a wide range in [Fe/H]. Stars in the 
family plotted in red (i.e., those with [a/Fe] < 0.2, see Fig. [2} are 
those traditionally assigned to the thin disk. Those in green are 
a-enhanced stars ([a/Fe] > 0.2) usually associated with the thick 
disk. Counterrotating stars unlikely to belong to either the thin 
or thick disks are shown in cyan. 



3.1 A chemical definition of the thin disk 

Fig. [T]suggests that stars with measured Fe, a and Eu abun- 
dances separate into two "families", one with high [Fe/H] 
and low [(a+Eu)/Fe] (shown in red) and the other with 
lower [Fe/H] and high [(a+Eu)/Fe] (green). This separation 
is confirmed in the top panel of Fig. [5] which is identical to 
Fig. Q] but for the index normalized to Fe. The [(a + Eu)/Fe] 
histogram, which should be relatively free from selection bi- 
ases, hints strongly at the presence of two distinct compo- 
nents. This separation highlights the well-known a enhance- 
ment of thick disk stars relative to the thin disk, accentu- 
ated by the addition of Eu to the a index. The data suggest 
that a plausible boundary between the two families may be 
drawn at [(a+Eu)/Fe]= 0.2, the value that roughly marks 
the "valley" in the histogram. 

Guided by this, we have chosen arbitrary but plausi- 
ble boundaries to define the two components in Fig. [1] and 
the top panel of Fig. [2] which we color in green and red, 
respectively. The simple criteria: (i) [Fe/H]> —0.7, and (ii) 
[(a+Eu)/Fe]< 0.2, distinguishes well one of the two families 
of stars. This family (coloured red) contains the Sun (which 
would be at the origin of the plot) and contains mostly stars 
associated with the thin disk as usually conceived. The mini- 
mum [Fe/H] boundary is suggested by the fact that the most 
metal-rich counterrotating star in the sample (certainly not 
a member of the thin disk) has [Fe/H]~ —0.7. 

The distinction between components blurs when includ- 
ing stars for which Eu abundances are not available (see 



Figure 2. Top panel: Same as Fig. [TJ but for the a+Eu index 
normalized to Fe. The presence of the two families of stars hinted 
at in Fig. [l] is confirmed by the [(a + Eu)/Fc] histogram. We 
identify the bottom sequence, shown in red, with the thin disk, 
and the top sequence, in green, with the thick disk. The conditions 
used to define the thin disk are (i) [Fc/H] > —0.7 and (ii) [(a + 
Eu)/Fe] < 0.2. Stars in the region defined by [Fe/H]> —1.5 and 
[(a + Eu)/Fc] > 0.2 are those traditionally associated with the 
thick disk. Note that this condition allows in some counterrotating 
stars (V < 0, shown in cyan) unlikely to be true members of the 
thick disk; the remainder are shown in black. Bottom panel: Same 
as top panel but for the index without Eu. Stars with available Eu 
are shown here using the same red and green colours as in the top 
panel. Stars in cyan are counterrotating (V < 0) stars, unlikely to 
belong to either the thin or thick disks. Note that counterrotating 
stars in the range —1.5 <[Fe/H]< —0.7 have lower [a/Fe] values 
than those of the thick disk. These stars appear to belong to a 
component ("D") separate from the thin and thick disk or the 
traditional stellar halo, which we link to tidal debris. See text for 
further details. 

bottom panel of Fig. [5J), but the sample more than doubles. 
In order to take advantage of this larger sample, we define 
the thin disk by the same criteria as above, namely 

• (i) [Fe/H]> -0.7; and 

• (ii) [a/Fe]< 0.2. 

We emphasize that the criteria above are purely chemi- 
cal. This differs from the traditional practice of selecting 
thin disk stars by their kinematics and might therefore in- 
clude stars not expected in kinematic definitions of the "thin 
disk" . We shall hereafter call them "thin disk stars" , but this 
caveat should be kept in mind when comparing our results 
with other work on the topic. 

Both panels of Fig. [2] show the wide range in [Fc/H] 
spanned by thin disk stars, as well as the strong correlation 
between the abundance ratio and metallicity in this popula- 
tion; [a/Fe] and [(a+Eu)/Fe] decrease slightly but steadily 
with increasing [Fe/H]. 

This trend suggests two possible interpretations. In 
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Figure 3. Mean rotation velocity and velocity dispersion as a 
function of [Fe/H] for stars identified as belonging to the thin 
disk component in Fig. [2] Filled circles correspond to stars with 
measured Eu abundances (top panel of Fig. open circles to 
all stars in the region labelled "Thin" in the bottom panel of 
Fig. [2] Error bars indicate one-sigma bootstrap error estimates. 
Note that, although they span more than a decade in [Fe/H], the 
average rotation velocity and velocity dispersion of thin disk stars 
is roughly independent of metallicity. The open squares shows 
the <r-[Fe/H] correlation for all stars in our sample and illustrates 
the familiar increase of the total velocity dispersion ( <r) with 
decreasing [Fe/H]. 



standard chemical evolution models it is reminiscent of a 
self-enriched population whose star formation timescale is 
long compared with the lifetime of stars that end their 
lives as supernovae type la (SNIa). These supernovae return 
mostly Fe and little a to the ISM, and therefore the ratio 
[a/Fe] of successive generations of stars declines steadily as 
[Fe/H] increases. 

In this interpretation, the formation of the metal-poor 
tail of the thin disk precedes that of the metal-rich one 
since [Fe/H] is assumed to march roughly monotonically 
with time. One difficulty with this interpretation is that 
the wide range in metallicity spanned by the thin disk 
(-0.7 <[Fe/H]< +0.4; see top panel of Fig. [2} requires a 
protracted star formation history and implies a strong cor- 
relation between age and metallicity. This seems at odds 
with the relatively weak age-[Fe/H] relation observed for 
thin disk stars. One must note here that stellar ages are no- 
toriously difficult to estimate, and that the strength of the 
predicted correlation depends on assumptions about infall 
and star formation efficiency that are difficult to pin down 
accurately. Despite these caveats, a number of authors have 
concluded that the lack of a strong age-metallicity relation 
signals the need to consider a different mechanism to exp lain 
the properties of the thin disk (see, e.g., lHavwoodll2008l . and 
references therein). 

These difficulties are circumvented in a second, differ- 



ent interpret ation of the properties of th e local thin disk. 
As shown bv lSchonrich fc Binnevl (|2009bh . radial migration 
can populate the vicinity of the Sun with "thin disk" stars 
that formed either inside or outside the Solar circle. The 
properties of local thin disk stars might not trace the chem- 
ical history of the solar neighbourhood, but rather reflect 
the spread in birth radii of such stars. Given the metallicity 
gradient in the Galaxy, this could explain naturally the wide 
range in metallicity of thin-disk stars found today in the so- 
lar neighbourhood. The same modelling predicts that, once 
a steady state is reached, additional enrichment at different 
radii would be balanced by the infall of fresh metal-poor gas 
from the IGM. Because of higher densities, star formation 
progresses nearer the center more intensely relative to the 
amount of available gas, leading to lower metallicities and 
higher [a/Fe] ratios with increasing Galactocentric birth ra- 
dius. These gradients are today reflected locally in the ob- 
served correlation between [Fe/H] and [a/Fe]. 

One disadvantage of this scenario is that it requires a 
steeper met allicity gradient in the ISM than is usually ac- 
cepted (see iMaciel et af] 120071 . and references therein), as 
well as careful balancing of the timescales of star formation, 
radial migration, and gas infall at different radii in order to 
explain the observed slope of the [a/Fe]- [Fe/H] correlation. 
Another potential challenge may arise from the fact that at 
least some metal-rich stars in th e inner Galaxy seem to be 
a-enriched (jMelendez et al.ll200Sf ) whereas those in the solar 
neighbourhood are not. 

One potential test of the migration scenario was dis- 
cussed in the Introduction, and concerns the prediction 
of a negative correlation between metallicity and rota- 
tion speed in the thin dis k. Reading off from Fig. 4 in 
ISchonrich fc Binnevl (|2009bl ), the mean rotation speed is ex- 
pected to change from roughly 250 to 160 km/s as the metal- 
licity in the thin disk increases from [Fe/H]= —0.7 to +0.5. 

As shown in Fig. [3] this prediction finds some support 
in the data. Filled circles here refer to the a+Eu data shown 
in the top panel of Fig. [5] whereas open symbols correspond 
to the much larger sample of stars without Eu data (bottom 
panel of Fig. [2]). The mean rotation speed of the thin disk 
(defined by the above conditions) decreases somewhat with 
[Fe/H] in the a+Eu dataset but th e trend is weaker than 
expected from ISchonrich fc Binnevl l |2009bl ) and the error 
bars are large. 

Furthermore, when the larger sample without Eu is con- 
sidered (open circles) the mean rotation speed is actually 
independent of metallicity. It appears as if the mean rota- 
tion speed, as well as the velocity dispersion, remain roughly 
constant over the whole range in [Fe/H] spanned by the thin 
disk. The constant velocity dispersion of the thin disk should 
be cont rasted with the fa miliar increase in a with decreasing 
[Fe/H] (|Stromgrenll 19871 ); this is shown by the open squares 
in Fig. [3] which corresponds to all stars in our compilation. 

The kinematic invariance of the thin disk with metallic- 
ity is a somewhat surprising result for either of the two sce- 
narios advanced above. The absence of a (V)-[Fe/H] correla- 
tion disfavours the migration scenario; on the other hand, in- 
situ self-enrichment scenarios would predict the velocity dis- 
persion to rise with decreasing [Fe/H] since metal-poor stars 
would be, on average, older, and therefore exposed for longer 
to gravitational heating by inhomogeneities in the Galactic 
potential. The data in Fig.[3]disfavours such gradual heating; 
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Figure 4. Top panel: The open histogram shows the distribu- 
tion of the rotation velocity (V component) for all stars in the 
sample. (Units are number of stars per bin.) The shaded his- 
togram corresponds to those in the thin disk, as defined in the 
bottom panel of Fig. [2] The remainder are shown by the open 
histogram. Bottom panel: V-distribution of stars in the metal- 
poor tail, i.e., [Fe/H]< —1.5. This illustrates the velocity distri- 
bution of the slowly-rotating, dynamically-hot stellar halo. The 
blue curve illustrates the best fitting gaussian, with (V) = —60 
km/s and ay = 144 km/s. Middle panel: V-distribution of stars 
with [Fc/H]> —1.5 but excluding the thin disk. The distribution 
is strongly non-gaussian, and shows two well-defined peaks; one 
at V ~ km/s and another at V ~ 160 km/s. The latter is well 
traced by stars in the "thick disk" component identified in the 
top panel of Fig. [2] (with Eu; filled green histogram). Including 
stars without Eu measurements but of comparable [a/Fc] ratios 
(the "Thick" region in the bottom panel of Fig. [2j results in the 
shaded green histogram. Note that stars near the V ~ peak 
are almost exclusively those in the debris ( "D" ) region of Fig. [2] 
They define a non-rotating, dynamically-cold component distinct 
from either the thick disk and the stellar halo. Both the debris 
and thick disk components arc relatively cold dynamically; the 
V-distribution can be well approximated by the sum of two gaus- 
sians with similar velocity dispersion, of order ~ 40 km/s (see 
green and magenta solid curves). 



if present, the heating mechanism must operate promptly 
and saturate quickly, as has be e n suggested in the past 
llStromgrenlll987l; iFreemanl Il99ll ; lOuillen fc Garnettl l200d : 
ISoubiran et alJfeoOct l. Inefficient heating would actually be 
easier to reconcile with theoretical models, which have strug- 
gled to explain the rapid increase in velocity dis persion 
with age inferred from earlier observation s (see, e.g . , IWielen 



2009) 



19771; iLacevll 19841 ; I Jenkins fc Binnevlll990l ; lAumer fc Binnev 



3.2 A chemical definition of the thick disk 

If our interpretation is correct, then the increase in velocity 
dispersion with decreasing metallicity for stars in the vicin- 



ity of the Sun must result from the increased prevalence of 
the thick disk at low metallicity. Indeed, as mentioned in 
Sec. [T] it is generally agreed that the thick disk is metal- 
poor; lags the thin disk in rotation speed; and has a higher 
velocity dispersion. This is shown in the top panel of Fig. 21 
where we compare the distribution of the rotation speed (V 
component) of all stars in our sample with that of the thin 
disk. 

The V distribution of stars not in the thin disk is com- 
plex, and hints at the presence of distinct dynamical compo- 
nents. It shows two well defined peaks, one at V~ 160 km/s 
and another at V ~ km/s, as well as a tail of fast coun- 
terrotating stars at highly-negative values of V. The first 
peak corresponds to a rotationally-supported structure: the 
traditional thick disk. Stars belonging to this rotating com- 
ponent seem to disappear from our sample when only metal- 
poor stars with [Fe/H] < —1.5 are considered (bottom panel 
of Fig. 2]). The latter trace the classical, kinematically-hot, 
metal-poor stellar halo: the V distribution is consistent with 
a gaussian with velocity dispersion av ~ 144 km/s (shown 
in the bottom panel with a blue curve). 

Interestingly, the V distribution of non-thin-disk stars 
with [Fe/H] > —1.5 (middle panel of Fig. |4| shows even 
more clearly the double-peak structure noted above. The 
peak at V ~ 160 km/s is well traced by the stars identified 
with the thick disk in the a+Eu panel of Fig. [2] shown by 
the solid-shaded histogram in the middle panel of Fig. 0] 

Indeed, the peak is traced almost exclusively by stars 
with high values of [a/Fe]. This is shown by the shaded 
green histogram, which corresponds to all stars in the re- 
gion labelled "Thick" in the bottom panel of Fig. d The V 
distribution of these stars is well approximated by a gaus- 
sian with (V) = 145 km/s and av = 40 km/s that accounts 
for nearly all stars with V > 100 km/s. 

The association between the thick disk and "high-a" 
stars is reinforced by inspecting the location of counterro- 
tating stars in Fig. [2] (shown in cyan). These stars, which 
clearly do not belong to a rotationally-supported structure 
like the thick (or thin) disk, are evenly distributed among 
stars with [Fe/H] < —1.5 but shun the "high-a" region in 
the range —1.5 < [Fe/H] < —0.7. Of the 38 counterrotating 
stars in our sample with —1.5 < [Fe/H] < —0.7, only 3 lie 
above the dotted line that delineates the "Thick" region in 
the bottom panel of Fig [2] 

It is tempting therefore to adopt a purely chemical def- 
inition of the thick disk in terms of [Fe/H] and [a/Fe]: 



-0.7; 



(i) [Fe/H] > -1.5; 

(ii) [a/Fe] > 0.2- ([Fe/H] +0.7) /4 for -1.5 < [Fe/H] < 
7; 

(iii) [a/Fe] > 0.2 for [Fe/H] > -0.7. 



If ou analysis is correct, then the thick disk emerges 
as a chemically and kinematically coherent component that 
spans a wide range in metallicity ( — 1.5 <[Fe/H]< —0.3) and 
contains mainly stars highly enriched in a elements. Stars 
in this component have ((U),(V),(W»= (10, 133, -2) km/s, 
and (<tu ,o~v ,o~w)= (95,61,61) km/s. 

The average V lags below the ~ 160 km/s peak because 
the chemical definition allows in a few stars with discrepant 
velocities (some with negative V) more closely associated 
with the V ~ peak than with a rotationally supported 
structure like a thick disk. This suggests that the chemical 
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definition of the thick disk proposed above is not perfect, 
and that it includes spurious stars belonging to a different, 
non-rotating component. We turn our attention to that com- 
ponent next. 



3.3 Tidal debris in the stellar halo 

The last feature of note in the top and middle panels of 
Fig. [4] is the peak centered at V ~ 0. These are stars with 
no net sense of rotation around the Galaxy and with a sur- 
prisingly low velocity dispersion. (The gaussian fit to the 
low-V tail shown in magenta in the middle panel of Fig. [4] 
has a dispersion of just 40 km/s.) Their vertical (W) velocity 
distribution (computed for stars with V < 50 km/s to better 
isolate the pealQ) is also indistinguishable from that of the 
thick disk component discussed above, with (W) ~ km/s 
and aw ~ 72 km/s. On the other hand, their U velocity 
dispersion is large (~ 180 km/s). Further, the U distribu- 
tion (not shown) is rather flat, with hints of two, symmetric 
peaks, one at U ~ —250 km/s and another at U ~ +250 
km/s. 

As discussed bv lMeza et all (|2005l ). this is the kinematic 
structure expected for a tidal stream originating in a dwarf 
galaxy whose orbital plane at the time of disruption was 
coincident with the plane of the Galaxy. The low V and 
W velocity dispersions are easily explained in this scenario, 
since such stream would be, locally, kinematically cold and 
confined to the orbital plane. The double-peaked U distri- 
bution also arises naturally in this scenario if stream stars 
have apocentric radii outside the solar circle and pericen- 
tric radii inside the solar circle. Stream stars in the solar 
neighbourhood are therefore either going to their apocenter 
with large, positive U or coming from their apocenter, with 
symmetric — U. 

This suggests that most non-thin-disk stars with 
[Fe/H]> —1.5 and low, but still enhanced relative to solar, 
[a/Fe] (previously thought to belong to the classical halo) 
belong to this new component. Guided bv lNissen fc Schuster! 
(|2010r i. we inspect the Na and Ni content of such stars for 
supporting evidence of this conclusion. This is shown in 
Fig. [5] where we show [Na/Fe] vs [Ni/Fe] for all stars in 
our sample with [Fe/H] < —0.7. Blue squares correspond to 
the very metal poor stars in our sample ([Fe/H] < —1.5); 
green open circles are stars in the a-rich "Thick" region of 
Fig. [2] and magenta filled circles denote stars in the a-poor 
"debris" ( "D" ) region of Fig. The three groups separate 
clearly in the Na-Ni plane, supporting our claim that the 
"debris" component is truly distinct from the thick disk and 
from the metal-poor "classical" halo. 

Using the same [Fe/H] and [a/Fe] parameters as above, 
we can characterize "debris" stars in the [a/Fe] vs [Fe/H] 
plane (region labelled "D" in Fig. O by 

• (i) -1.5 < [Fe/H] < -0.7; and 

• (ii) [a/Fe] < 0.2 - ([Fe/H] + 0.7) /4. 

Our conclusion agrees with that of iNissen fc Schuster] 
(2010). who studied a large spectroscopic sample of "halo" 
stars and argued, in agreement with our analysis, that most 



1 Lifting this restriction has only a small impact on the disper- 
sion; aw ~ 90 km/s for all stars in the "D" region of Fig. [2] 
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Figure 5. [Na/Fe] vs [Ni/Fe] correlation for all stars in our sam- 
ple with [Fc/H]< —0.7. Blue squares correspond to [Fc/H]< —1.5 
"classical halo" stars; green open circles to a-rich "thick disk" 
stars with —1.5 < [Fe/H] < —0.7, and magenta filled circles to 
stars in the "debris" ("D") region of Fig. [2] Note how these three 
components separate neatly in the Na vs Ni plane and, in partic- 
ular, the tight scatter around the mean trend of "debris" stars. 
This strengthens our conclusion that they correspond to three 
components of distinct origin. 



metal-rich "a-poor" halo stars are indeed tidal debris from 
disrupted dwarfs. Our debr is population is remin iscent of 
the population identified bv lMorrison et all (120091 ). 



4 SUMMARY AND DISCUSSION 

The preceding analysis suggests that apportioning the var- 
ious components of the Galaxy according to purely chem- 
ical criteria is both possible and fruitful. The definition of 
the thin disk in the [(a+Eu)/Fe] vs [Fe/H] plane is par- 
ticularly straightforward, and suggests that the kinematics 
of the thin disk is invariant with metallicity. This is an in- 
triguing result unexpected in migration-based scenarios for 
the chemo-dynamical evolution of the thin disk. It implies 
that the familiar in crease in velocity dispersion with de- 
creasing metallicity jStromgrenl 1 19871 ) is the result of the 
increased prevalence of the thick disk at lower metallicities, 
rather than of the sustained operation of a dynamical heat- 
ing mechanism. If confirmed, the kinematic invariance of the 
thin disk with metallicity will place strong constraints on the 
formation of the Galactic disk and on the role of accretion 
events, in situ formation, and/or migration. 

The "thick disk" can also be charted in the [a/Fe] vs 
[Fe/H] plane. As reported in earlier work, it seems to contain 
mainly stars highly enriched in a elements. It shows as a 
separate dynamical component in rotation speed, with the 
bulk of its stars rotating at V ~ 160 km/s. A simple criterion 
in a content isolates most of these stars, although a few 
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outliers with nearly zero, or negative, V velocities are also 
included. The latter might very well be contaminants from 
a different population that a crude boundary in the a-Fe 
plane is unable to weed out. 

The inclusion of additional heavy elements in the defin- 
ing criteria might enable a cleaner characterization of the 
thick disk. If that were possible, questions such as whether 
the thick disk shows evidence of self-enrichment, or whether 
correlations between metal content and velocity dispersion 
are present, could be addressed. This would allow us to 
distinguish between migration and accretion models and, 
among the latter, between those where the bulk of thick 
disk stars were either accreted or simply stirred. 

A substantial fraction of stars in the range —1.5 < 
[Fe/H] < —0.7 seem to belong to a dynamically-cold, non- 
rotating component with properties consistent with those 
of a tidal stream. These are mainly stars of low-a content, 
comparable in that regard to individual stars in many of the 
satell i te companions of the Milky Way (see, e.g.. lVenn et al.l 
|2004| ; iTolstov et all l2009h : the low- [Fe/H], a-poor region 
should be a good hunting ground for the remnants of ac- 
cretion events. 

The V-distribution associated with the stream or "de- 
bris" seems to peak at slightly negative V (~ —10 km/s). 
This implies that most stars in the stream have angu- 
lar momenta similar to that of the globular cluster uCen 
jDinescu et all 1 19991 ). It is therefore tempting to associate 
these stars with the parent galaxy of this massive clus- 
ter, long suspected to be the survivor of t he disruption o f 
a dwarf galaxy in the Galactic potential |Freemanl 1 19931 ). 
Supporting evidence comes from cjCen's low vertical ve- 
locity, as well as from the overlap in m etallicity between 
q>Cen and stars in the pu tative stream (see lMeza et al ] |2005l ; 
iNissen fc Schusterll2010l . for a full discussion). 

Although the orbital parameters of the stream might 
suggest a link with the globular cluster LuCen, we caution 
that evidence of a true relation between those stars and 
LjCen is circumstantial at best. For example, t he metallicity 
distribution in wCen peaks at [Fe/H] = -1.7 dSmithl 120041) 
and have distinct elemental abunda nces (|Norris fc Da Costal 
1 19951 ; IJohnson fc Pilachowskill2010l ) , so it might be useful to 
seek evidence for a stream in stars of similar metallicity and 
peculiar abun dance ratios when larger sam ples become avail- 
able (see, e.g.. IWvlie-de Boer et al" |2010| . for some progress 
in this direction). 

One corollary of this finding is that very few stars with 
[Fe/H] > —1.5 in our sample seem to belong to the classi- 
cal, dynamically-hot stellar halo. Recent work has suggested 
that halo stars in this "metal-rich" tail might a ctually belong 
to a d istinct "inner halo" component (see, e.g.. lCarollo et ah! 
[20071 . and references therein). The connection between this 
and the stream we advocate above is, however, unclear, not 
least because the "inner halo" component is reported to co- 
rotate with the Sun, whereas our putative stream counter- 
rotates slowly around the Galaxy. 

The "inner halo" identification typically relies mainly on 
estimates of [Fe/H] and lacks information on the abundance 
of individual elements. On the other hand, our sample is 
relatively small in comparison and has potentially a num- 
ber of selection biases. Thus the possibility that our stream 
constitutes a small subset of the inner halo remains. 

It might be possible to test these ideas by examining 



other abundance ratios, such as the neutron-capture ele- 
ments [Ba/Y] or [Ba/Eu]. This is because these elements 
are contributed in different amounts by AGB stars that 
undergo slow neutron capture nucleosynthesis during the 
thermal pulsing stages and by massive stars that undergo 
rapid neutron capture nucleosynthesis during supernova ex- 
plosion. Indeed, variations in those ratios have already been 
seen in individual stars of dwarf galaxies (see Figure 14 in 
ITolstov et"al"1l2009l ). 

Overall, our success in dividing and assigning stars of 
the solar neighbourhood to families of distinct chemistry 
and kinematics seems to favour models where accretion 
events have played a significant role in the formation of the 
Galaxy rather than models, such as those based on migra- 
tion, where secular evolutionary mechanisms rule. We has- 
ten to add, however, that the criteria to separate compo- 
nents proposed here are imperfect, and that our conclusions 
are based on small and heterogeneous samples. These sam- 
ples likely conceal a number of biases which can only be 
revealed and lifted by a concerted effort to survey a large, 
volume-limited, kinematically-unbiased sample of stars with 
the high-resolution spectra needed to measure the abun- 
dance of individual heavy elements. The planned HERMES 
survejQ should be a major first step in this regard. This 
undoubtedly ambitious endeavour would allow us to dissect 
chemically and kinematically the solar neighbourhood and 
to learn the true provenance of the many stellar families that 
today call this small place of the Galaxy home. 
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